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We explore the non-equilibrium response of Chern insulators. Focusing on the Haldane model, we 
study the dynamics induced by quantum quenches between topological and non-topological phases. 

A notable feature is that the Chern number, calculated for an infinite system, is unchanged under the 
dynamics following such a quench. However, in finite geometries, the initial and final Hamiltonians 
are distinguished by the presence or absence of edge modes. We study the edge excitations and 
describe their impact on the experiment ally-observable edge currents and magnetization. We show 
that, following a quantum quench, the edge currents relax towards new equilibrium values, and that 
there is light-cone spreading of the currents into the interior of the sample. 
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Topological phases of matter display many striking fea¬ 
tures, ranging from the precise quantization of macro¬ 
scopic properties, to the emergence of fractional excita¬ 
tions and gapless edge states. An important class of topo¬ 
logical systems is provided by the so-called Chern insu¬ 
lators realized in two-dimensional settings [1]. A famous 
example is the Haldane model [2], which describes spin¬ 
less fermions hopping on a honeycomb lattice. The Hal¬ 
dane model exhibits both topological and non-topological 
phases, and its behavior is closely related to the integer 
quantum Hall effect. Recent advances using ultra cold 
atoms [3-7] have led to the experimental realization of 
the Haldane model [8]. Proposals also exist for realizing 
other states of topological matter using cold atoms [9]. 

A fundamental characteristic of topological systems 
is their robustness to local perturbations, making them 
ideal candidates for applications in metrology and quan¬ 
tum computation. However, much less is known about 
their dynamical response to global perturbations and 
time-dependent driving. This issue is of relevance in a 
variety of contexts, ranging from the time-evolution and 
controlled manipulation of prepared topological states, 
to the dynamics of topological systems coupled to their 
environment. Understanding the impact of topology on 
the out of equilibrium response is crucial for further de¬ 
velopments, and is the motivation for this present work. 
For recent progress in this direction see Refs. [10-17]. 

In this manuscript we investigate the non-equilibrium 
dynamics of the paradigmatic Haldane model. In partic¬ 
ular, we consider quantum quenches and sweeps between 
topological and non-topological phases. Key questions 
that we will address include: What happens to the topo¬ 
logical properties on transiting between different phases? 
What happens to the edge excitations following a quan¬ 
tum quench? How do the topological characteristics in¬ 
fluence the non-equilibrium dynamics? 

Model .— The Haldane model describes spinless 
fermions hopping on a honeycomb lattice with both near¬ 
est and next nearest neighbor hopping parameters. The 



Figure 1. (a) Phase diagram of the Haldane model obtained 
from the low-energy Dirac fermion representation, showing 
topological ( v — ±1) and non-topological phases (v — 0) [2]. 
We consider quantum quenches and sweeps between different 
regions of the phase diagram, as illustrated by the arrows, (b) 
The low-energy spectrum of the Haldane model is described 
by excitations around two Dirac points. After a quench, car¬ 
riers in the lower band are excited to the upper band. 

Hamiltonian is given by [2] 

H = hJ2 (ctcj.+hx.) 

(i,j) 

+ t 2 ^ [e^^clc- + h.c.^ 

((id)) 

+ M^hi - M^hi, ( 1 ) 
ieA ieB 

where the fermionic operators obey the anticommutation 
relations = 5ij and hi = cjdi. Here, (i,j) and 

((i,j)) indicate the summation over the nearest and next 
to nearest neighbor sites respectively, and A and B la¬ 
bel the two sub-lattices. The phase factor cpij = =b ip 
is introduced in order to break time-reversal symmetry 
and is positive for anticlockwise next to nearest neighbor 
hopping. The energy off-set ±M breaks spatial inversion 
symmetry. The phase diagram of the Haldane model is 
shown in Fig. 1 (a); following Ref. [2] we assume that 
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|^2All < 1/3 so that the bands may touch, but not over¬ 
lap. 

For ^ 2 ,M < £i, the Hamiltonian (1) has a linear dis¬ 
persion near the six corners of the hexagonal Brillouin 
zone, but only two of these are inequivalent. As a result, 
close to half-filling, the low-energy description is given by 
the sum of two Dirac Hamiltonians 


/ m a c 2 — cke ia0 \ 

\—cke ~ ia0 —m a c 2 J ’ 


( 2 ) 


where a = ±1 label the Dirac points. Here, c = 3t\/2h 
is the effective speed of light, /cexp(i$) parameterizes the 
2D momentum (k x ,k y ) and m a = (M — 3\/3<^£2 smp)/c 2 
is the effective mass [2]. The topological phases have a 
non-vanishing Chern number v [1, 2, 18, 19]. For a state 
1^) this is defined by the integral of the Berry curvature 
over the 2D Brillouin zone 



where ft = d kx A ky - d ky A kx and A ky = i (ip | d ky \ ip) is 
the Berry connection. For the ground state of the Hal¬ 
dane model v G ±1,0. This may be decomposed into 
contributions from the two Dirac points as v = z/ + + v _, 
where v a = — ^sign(m a ) G ±1/2. The boundaries of 
the topological phases correspond to the locations where 
m± changes sign. They are thus given by Mft^ = 
±3\/3sin<£, and are independent of t\\ see Fig. 1. 

Quantum Quenches .— In order to gain insight into the 
non-equilibrium dynamics of the Haldane model, we con¬ 
sider quantum quenches between different points (M, p) 
on the phase diagram shown in Fig. 1, for fixed values of 
t\ and £ 2 • At time t = 0, we prepare our system in the 
ground state with parameters (M 0 ,(^o)- At half-filling 
our initial state fills the lower band. We then abruptly 
change the parameters of H to (M, </?), and allow the 
system to evolve unitarily under the action of this new 
Hamiltonian. In general, this will lead to a non-trivial 
occupation of both the lower and the upper bands. 

We begin by examining the non-equilibrium response 
of the effective Dirac Hamiltonian H = H+ + H _. Since 
v = [sign(ra + ) — sign(ra_)], quenching between dif¬ 
ferent phases corresponds to changing the sign of one 
or both of the masses m a . For a given Dirac point, 
such changes will lead to a re-distribution of carriers be¬ 
tween the two bands. For a 0 independent superposition, 
\ip a (k)) = \l a {k)) +b a {k)e~ iE ^ t \u a (k)), 

the Chern number is formally given by 


Va{t) = —asign m a Q _ |& a (0)| 2 ^ 

- |6 a (oo)||a a ( 00 )| cos[(.E“(oo) - E l a (oo))t + 5], (4) 

Here, a a (k) and b a (k) are complex c- numbers, S = 
arg(a a (oo)) — arg(6 a (oo)), and E l ^ u (k) are the energies 



Figure 2. Probability of occupying the upper band for a sin¬ 
gle Dirac point (a = —1) with c = 1, following a quench of 
m a . Sign-preserving quench, m_ = — 1 —»> m'_ — —0.1 (solid 
line) and a sign-changing quench, ra_ = — 1 —> m'_ = 0.1 
(dashed line). In both cases, 6_(oo) = 0, corresponding to the 
time-independence of i/_(£) using Eq. (4). The sign-changing 
quench yields |6_(0)| 2 » 1, but the contribution to v~(t) in 
Eq. (4) is compensated by the change in sign of m_. As a 
result, v- is unchanged from its initial value. 


in the lower and upper bands. In general, u a (t) is time- 
dependent, and differs from its ground state values ±1/2. 
However, the time-dependence only enters via the super¬ 
position coefficients evaluated at k = oo. An explicit 
computation shows that b a { oo) = 0, following a quan¬ 
tum quench; see Fig. 2 and the Supplemental Material. 
In addition, b a ( 0) = 0, ±1, so the potential modifica¬ 
tion of v a is compensated by the change in sign of m a . 
As a result the Chern number is unchanged from its ini¬ 
tial value , even if one quenches between different phases. 
Similar results may also be obtained for a linear sweep, 
m a(t) = £/t; see Supplemental Material. 

Preservation of Chern Number. — An intuitive way 
to understand the persistence of v following a quantum 
quench is in terms of spin-textures in momentum space. 
The Dirac Hamiltonian in Eq. (2) can be recast as an 
effective spin in a k-dependent magnetic field, h Q! (k). 
Explicitly, H a { k) = — h a (k) • cr/2, where a are the 
Pauli matrices. In equilibrium, the topological phases 
with u a = ±1/2 correspond to meron spin configura¬ 
tions which wind on the upper (lower) half-sphere [20]. 
Following a quantum quench, the spins precess in the ef¬ 
fective magnetic field of the new Hamiltonian, preserving 
the topological characteristics of the initial spin configu¬ 
ration. A similar argument may also be applied to the 
Haldane model (1) in k-space. Indeed, one expects the 
preservation of topological invariants under time evolu¬ 
tion to be a general feature for non-interacting fermions 
in a periodic system, where each k-state evolves uni¬ 
tarily under some Hamiltonian H( k), provided H( k) is 
smoothly varying in k-space. 

Edge States. — In the above discussion we have demon¬ 
strated that the value of v is unchanged as one quenches 
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Figure 3. Energy spectrum of the Haldane model obtained 
by exact diagonalization on a finite-size strip of width N = 
20 unit cells with armchair edges. We take periodic (open) 
boundary conditions along (transverse to) the strip and set 
£i « 1, = §, and M — 1. (a) Equilibrium population of the 

energy levels in the non-topological phase with p = fr. (b) 
Re-population of the levels after a quench to the topological 
phase with p — corresponding to the solid arrow in Fig. 1 
(a). The size of the dots is proportional to the probability of 
finding a particle in the mode. Post-quench, the filling of the 
edge states and the bands is non-trivial. 


dinal current flowing along the strip in the x-direction 
at a definite transverse ^/-position is therefore given by 
J n = (Jn) = T,m S (Jmns)- In Fi S- 4 ( a ) we P lot this cur - 
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and sweeps between different phases. However, there is 
a fundamental distinction between the topological and 
non-topological phases, due to the presence or absence 
of edge states in a finite-size sample [21]. In quenching 
between phases of different topological character, these 
edge states will either appear or disappear, depending on 
the direction of the quench. This is confirmed in Fig. 3, 
which shows the re-construction and re-population of the 
energy levels following a quench from the non-topological 
phase to a topological phase. It is readily seen that the 
edge states emerge and are populated as a result of the 
quench, in spite of the fact that v remains equal to zero in 
the absence of boundaries. Conversely, a quench from a 
topological phase to the non-topological phase eliminates 
the edge states, whilst v remains pinned at unity. 

Edge Currents and Orbital Magnetization .— Having 
examined the re-population of the edge states we now 
consider physical observables that depend on these states, 
including the edge currents and the orbital magneti¬ 
zation. We first consider these quantities in equilib¬ 
rium, which already display interesting features. We 
define the local current flowing through the site i by 
= —^^2jSji(tijc\cj — h.c.), where tij is the hop¬ 
ping parameter of the Haldane model between sites i 
and j, Sji is the vector displacement of site i from j, 
and the sum is over the nearest and next nearest neigh¬ 
bors. The site indices may be decomposed into the triplet 
{m,n,s} labeling the x and y positions of the unit cell 
and the sublattice index s = A, B. The total longitu- 


Figure 4. Equilibrium properties of the Haldane model on 
a finite-size strip as used in Fig. 3. (a) Total longitudinal 

current Jf along the strip as a function of the transverse 
spatial index n G 1,..., 20, for M = 0 and p = tt/3. (b) Edge 
currents corresponding to Jf with n — 1 (solid) and n — 20 
(dashed) for M — 0. The edge currents exhibit 7r-periodicity 
in p and vanish when p = 7r/2. (c) Orbital magnetization 
A4 as a function of p for M = 0. (d) Intensity plot of A4 
where the dashed lines correspond to the boundaries of the 
topological phases. Numerically we observe that vanishes 
on the loci M = d= sin p (solid) within the topological phases. 
The loci are fits to the numerical data (triangles) where M, IS 
0. The magnetization also vanishes on the vertical lines p — 
7T/2,7T, 37t/2,. .. as follows from symmetry considerations. 

rent within the topological phase for M = 0 and p = n/3. 
The presence of the counter-propagating edge currents 
is readily seen. In Fig. 4(b) we show the dependence 
of these edge currents on p. Somewhat surprisingly, the 
edge currents vanish within the topological phase, in spite 
of the presence of edge states in the spectrum. The edge 
currents are composed of counter-propagating contribu¬ 
tions which cancel at p = 7r/2; see Supplemental Ma¬ 
terial. Moreover, the longitudinal currents Jif exhibit 
7r-periodicity in p. This is a consequence of being at 
half-filling and occurs in spite of the fact that the Hamil¬ 
tonian and the current operator have a periodicity of 27r. 
To prove the 7r-periodicity in p we first note that both the 
Hamiltonian and the current operator change sign under 
the transformation M —M, p -T <^+7r, c mn A —>■ c mn ^, 
CmnB ~CmnB , thereby interchanging the upper and 

the lower bands. At half-filling, we fill only the lower 
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band, and it follows that J%(M,ip) = J^(—M, p + 7r). 
In addition, the current changes sign under the parity 
transformation x —»• —x. This interchanges the sublat¬ 
tices and corresponds to M 4 — M and ip —>* — <^. It 
follows that JZ(M,<p) = —cp) = J£(-Af», 

where in the last step we use the transformation proper¬ 
ties under time-reversal. Combining these relations, one 
obtains the 7r-periodicity in tp and the vanishing of the 
longitudinal currents for p = tt/2. 

Similar arguments also apply to the (lattice discretiza¬ 
tion of the) orbital magnetization: 

M = T/d 2 rr x (J(r)), (5) 

where J (r) is the local current density operator and A is 
the area. As shown in Fig. 4(c) this also vanishes within 
the topological phases and has 7r-periodicity in cp [22, 23]. 
Our numerical computations also reveal that the mag¬ 
netization vanishes on a sinusoidal locus M = =b sm(ip) 
within the topological phases; see Fig. 4(d). In addition, 
A4(M, (p) has extrema at M = 0 and on the topological 
phase boundaries, M = ±v / 3sin((^), for fixed (p . Away 
from half-filling, the particle-hole symmetry is broken 
and the periodicity of the currents and the magnetiza¬ 
tion is restored to 27r. The increase or decrease of the 
edge currents depends on the sign of the doping and the 
Chern index; see Supplemental Material. 

Dynamics of the Edge Currents. — Having discussed 
the equilibrium properties of the edge currents we now 
consider their response to quantum quenches. In Fig. 5 
we show quenches from the topological to the non- 
topological phase. The edge currents decay towards new 
values that are found to be numerically close to the equi¬ 
librium values of the post-quench Hamiltonian. This is 
in spite of the fact that the system is left in an excited 
state under unitary evolution, and that v remains pinned 
to unity in the absence of boundaries. Quenches from 
the non-topological to topological phases exhibit simi¬ 
lar behavior; see Supplemental Material. Further insight 
into the non-equilibrium evolution may be gleaned from 
the time-evolution of the longitudinal currents across the 
two-dimensional system. As shown in Fig. 6, the damped 
oscillations of the edge currents is accompanied by the 
light-cone spreading of the currents into the interior of 
the sample. It would be interesting to observe this dy¬ 
namics in experiment, which is in principle possible if 
local imaging is available [24]. 

Conclusions. — In this manuscript we have explored 
the non-equilibrium dynamics of the Haldane model. We 
have demonstrated that the Chern number is preserved 
in both quenches and sweeps between different regions of 
the phase diagram. However, the edge states may be re¬ 
constructed and re-populated leading to changes in the 
accompanying edge currents. Predictions for experiment 
include the vanishing of the equilibrium edge currents in 



Figure 5. Dynamics of the edge current J^{t) for N = 30 
(circles) and N = 40 (crosses) following a quantum quench 
between the topological phase and the non-topological phase 
with t\ = 1, t 2 = | and fixed <p — ir/3. Quenches from 
M = 1.4 to M = 1.6 (main panel) and from M = 1.4 to 
M — 2.2 (inset) showing that the edge currents approach 
new equilibrium values. For the chosen parameters, these are 
very close to the ground state expectation values of Jf? in the 
final Hamiltonian, as indicated by the horizontal lines. 
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Figure 6. Dynamics of the currents \Jf (t) | following a quan¬ 
tum quench from the topological to the non-topological phase 
for the parameters used in the inset of Fig. 5. The damped 
oscillations of the edge currents are clearly visible, as is the 
light-cone spreading of the currents into the interior of the 
sample, where c = 3t\/2h = 3/2 is the effective speed of 
light. The waves propagating from the two edges meet at time 
t ~ ( N/2)V3/2c ~ 5.77, leading to resurgent oscillations in 
finite-size samples; see Supplemental Information. 

the topological phases, and the light-cone spreading of 
the currents following a quantum quench. There are a 
wide variety of directions for further research, including 
the dynamics of the conductivity and its relation to the 
Chern number, and the role of decoherence via coupling 
to the environment. 
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SUPPLEMENTAL MATERIAL 


Haldane Model .— For completeness, let us recall some 
details of the Haldane model [2]. A key feature is the 
presence of time-reversal symmetry breaking induced by 
the complex second-neighbor hopping. In Haldane’s orig¬ 
inal paper this corresponds to a staggered magnetic field, 
where the local flux within a hexagonal plaquette is non¬ 
zero, but the total flux through a plaquette vanishes. In 
the experimental realization of Ref. [8] this is achieved 
by circular modulation of the lattice position. In Fourier 
space, the Hamiltonian is given by [2] 


H( k) = 2t 2 cos ip cos(k • b*) J I 

+ t\ ^^^[cos(k • a i)^ 1 + sin(k • a i)o 
M — 2t 2 sin ip ( sin(k • hi) 


+ 


a 3 , (6) 


where are Pauli matrices and I is the identity matrix. 
Here, ai,a 2 ,a 3 are the displacements from a B site to 
its three nearest-neighbor A sites. The b^’s are defined 
via cyclic permutations of bi := a 2 — a 3 and correspond 


to the second-neighbor displacements, ±b^. For M = 0 
and ip = 0, the two bands touch at the six corners of 
the Brillouin Zone. However, only two of these points 
are inequivalent, and we denote these by k a with a = d=. 
Expanding the Hamiltonian around these two corners, 
with U a = (n®, Ul) = k-k a , one obtains [2] 

H a = c(n* <7 2 - II 2 cr 1 ) + m a c 2 a 3 , (7) 

where m a and c defined following Eq. (2) and 

(U + in a) = ^X eik “' ai y^p- ( g ) 

In order to connect with the notations used in Eq. (2) we 
parameterize the complex number + in* as ke ia9 . 

Berry Phase. — In a gauge theory, invariant under a 
local transformation |?//(X)) := e~ lp |^(X)), parame¬ 
terized by X, one may define a covariant derivative 

|Z^(X)> = I^V’(X)) - IV’(X)) (V»(X) I 0^(X)>, (9) 

where |^(X)) (^(X)| projects out the parts of \d^{X)) 
that are not orthogonal to |^(X)). One may also intro¬ 
duce the Berry connection A M (X) = i (^(X) \d^(X.)) 7 
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the Berry phase </>, for a closed path T in parameter 
space e 1 ^ := e 1 $r dX ^ A »( x -\ and the Berry curvature 
= <A^( X ) I £W>(X)> - (D a tl>(X) I D^(X)). One 
may also generalize the Gauss-Bonnet theorem: 

v:=f = L[ dX^AdX ff ^ CT (X) ez, ( 10 ) 

Z 7 T Z 7 T Jm 2 

where M 2 is a closed orientable 2-manifold and v is the 
Chern number. 

Chern Number for Dirac Points .— For the Dirac 
Hamiltonian, given by Eq. (2) with a = — 1 say, the 
Berry connection is given by A jjT = i {ip \ dk'ip) and Af = 
i | V’) 5 where we parameterize the momentum-space 

2-manifold M 2 by (fc, 6). For a superposition of the form 


\ip(k, 0)) = a(k) \l(k, 9)) + b(k) | u(k, 9)), (11) 


where 1 1) and | u) correspond to the lower and upper band 
eigenstates, A jjT and Af are independent of 6. Asa result, 
the contribution of this Dirac point to the Chern number, 
v — V-\- H- V —, is 


1 f‘Z7T f‘CO f- OO 

v- = —\ d 9 d kSl~= dkdtAg = Ag\™ , 

Jo Jo Jo 

( 12 ) 

where D = dk A e — doA k . Using the explicit forms 

in _ (e~ i 9 f-(k,m-)\ , , _ / -e~ w f+(k,m-)\ 

f+(k,m-) )' W ~{ f-(k, m_) ) 

(13) 

where f ± (k,m) := (l ± yp=T 


, one obtains 


A e = \a(k)\ 2 f*(k,m-) + \b(k)\ 2 fl(k,m-) 

- j (a*(k)b(k) + h.c.). (14) 

2y k 2 + m?_ 

It follows that 


The probability of being in the upper band, |6_(fc)| 2 is 
plotted in Fig. 2. In particular, one obtains 6(oo) = 0 and 
6(0) = 0 (6(0) = ±1) for sign-preserving (sign-changing) 
mass quenches. The preservation of the contribution to 
the Chern number is discussed in the main text. 

Linear Sweep. — We extend the results for the Chern 
number preservation to linear, time-dependent, sweeps. 
We again consider a = — and set ra_ = ^ in Eq. (2), 
corresponding to a sign changing sweep over the interval 
t G [—oo, oo]. We prepare the system in the ground state 
at t = —oo and track the subsequent evolution [15]. At 
any instant of time the state can be written as a super¬ 
position over the eigenstates of H-(t): 

|^(fc, 0, t)) = a(k , t) |Z(fc, 6 , t )) + 6(fc, t) | u(k, 6 , t)}. (17) 
Using Eq. (12) and Eq. (14) one obtains 


v-(t) = sign t Q - \b{0,t)\ 2 ^ 





k 

+ 

lim — lim 


- (q*(fc, £)6(fc, t) + h.c.) 


k—Yoo k—> 0 


2 y k 2 + ^2 


(18) 

We now need to evaluate the coefficients a and 6 for k = 0 
and k = oo. For k = 0 the eigenstates of H-(t) are 


l*(0,t)> 

KM)) 


/ e “bK 1 - si s ni )\ 

\ bK 1 + si gni) / ’ 

f—e~ ie \J\0 + sign t)\ 

\ fiO- si s ni ) / 


(19) 

( 20 ) 


These are time-independent for t ^ 0, and for k = 0 
the system remains in the lower band up to t = CD. At 
t = 0 + the mass parameter changes sign, and at k = 0, 
overlaps completely with the upper band: 


A e (°°) = o I 1 “ (a*(oo)6(oo) + h.c.)] 


^V(O) 


16(0) | 2 m_ > 0 
|a(0)| 2 m_ < 0. 


(15) 

(16) 


Eq. (4) follows by reinstating the time dependence 
a(k) a(k)e~ lElt and b(k) b{k)e ~ lEut , for unitary 
evolution under the Dirac Hamiltonian. 

Quench .— We now consider a quench m_ —»• m'_. 
The system is initially prepared in the ground state of 
iJ_(m_), corresponding to the filled lower band, |/). Im¬ 
mediately after the quench we may decompose 1 1) into 
the eigenstates of the post-quench Hamiltonian: 


|<Z(0,0-)|u(0,0+))| = l. 

Henceforth, the k = 0 mode remains in the upper band. 
At k = oo the eigenstates are independent of time: 


|((oo, 



( 21 ) 

( 22 ) 


The system remains in the lower band for k = oo. Sum¬ 
marizing, one obtains 


a (k) = f-(k,m-)f-(k,m'_) + f + (k,m-)f+(k,m'_), 
b(k) = f+(k,m-)f-(k,m'_) - f-(k, m_ )/+(&, m'_). 


IKM)I = 


0, t < 0 
1, t > 0 


and 16(oo, t) | = 0 \/t. (23) 
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This parallels the situation for the quench protocol as 
shown by the dashed line in Fig. 2. The preservation of 
V- follows by substituting Eq. (23) into Eq. (18). 

Preservation of Chern Number. — Having established 
the preservation of v using the low-energy Dirac Hamil¬ 
tonian, we examine the non-equilibrium response of the 
Haldane model. By recasting Eq. (3) in the form 

v=fcf d (24) 

Z7T Jd BZ 

we may focus on the time-dependence of the Berry con¬ 
nection. In general, A k „ = i{ip\d k Jxf) + i(ip\d k ^) = 
(^|[^, H] |^), or equivalently, A kfM = {^\{d k ^H)\^). Ex¬ 
panding the initial state \^f) = u c 1 e~ lE A |y) in 

terms of the eigenstates |y) of the final Hamiltonian 


m 



A k ^{k x ,k y ) = ^2 c 7 Cy(7 , |a fe(iJ ff|7)e* (B V (25) 

7=/,n 


In general, this is time-dependent. However, using the 
symmetries of the final Hamiltonian, the components of 
A k[x along the Brillouin zone boundary occur in equal 
pairs and cancel in the line integral for z>. For example, 
within the upper and lower triangles depicted in Fig. 7, 
H(k x ,ky) = H*(k x , —ky). Periodicity in k- space ensures 
that H(k x ,k y ) = H(k x ,—k y ) along the corresponding 
zone boundaries, so that H is real on these segments. It 




Figure 7. First Brillouin zone of the Haldane model, (a) 
In each of the triangles the Berry connection and the Berry 
curvature are time-dependent. However, the Chern number 
is given by the line integral of the Berry connection along 
the zone boundary, as indicated by the arrows. Following a 
quench, the time-dependent contributions to v from opposite 
sides of the boundary cancel each other, (b) Time derivative 
of the Berry curvature for M = 1, t± 5 = 1 and 1 2 — 1/3 
following a quench from the non-topological phase with p — 
7 r/6 to the topological phase with ip = 7r/3. Although Q ^ 0, 
numerical integration over the Brillouin zone confirms that 
v = 0. 


follows that A kx (k x , 2tt/3) = A kx (k x , — 2tt/3) for k x G 
(— so these two contributions to z> cancel. 

Edge Currents and Orbital Magnetization. — In order 
to examine the behavior of the edge currents we consider 


Figure 8. The strip geometry used for our finite-size compu¬ 
tations. The strip is N cells wide with arm chair edges and 
has periodic (open) boundary conditions in the longitudinal 
(transverse) directions. Each unit cell is labeled by two in¬ 
dices m and n, and contains two sites belonging to the A or 
B sublattices. The phase p of the Haldane model is taken as 
positive for anticlockwise next to nearest neighbor hopping. 


the Haldane model on a finite-size strip with armchair 
edges and periodic (open) boundary conditions along 
(transverse to) the strip. The geometry we use is shown 
in Fig. 8. The numerical computations in the main text 
are performed on strips of width N = 20, 30 or 40 unit 
cells. We define the longitudinal current Jif as the to¬ 
tal current flowing along the n-th row of the strip where 
n G 1, ...,7V, as shown by the shaded area in Fig. 8. 
Exploiting the periodic boundary conditions along the 
strip, J* = J n (k x ). As shown in Fig. 9, the terms in 
this summation appear with both positive and negative 
signs. In particular, the edge currents, (Jf) and (Jjy-), re¬ 
ceive opposite contributions which perfectly cancel when 
(p = 7r/2. Doping the system with particles or holes 
breaks this symmetry and restores the 27r-periodicity, as 
shown in Fig. 10. The effect of doping is to change the 
relative contributions of the bulk and the edge states to 
the total edge currents; for particle (hole) doping the edge 
(bulk) states contribute more to the overall edge current. 

Zeros of the orbital magnetization also occur within 
the topological phases. As shown in Fig. 11, A d(M,p) 
has extrema at M = 0 and on the boundaries of the 
topological phases at M = Tv^sin^. Numerically we 
observe that the zeros of Ai occur on the sinusoidal loci 
M = dz sin p. Once again, the effect of doping is to 
change the relative contributions of the bulk and the edge 
states to the orbital magnetization; for particle (hole) 
doping the edge (bulk) states contribute more to Ai. 

Dynamics of the Edge Currents. — Following a quan¬ 
tum quench, we observe that the edge currents relax to 
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A 'X 

Figure 9. Momentum space contributions to the equilibrium 
currents along the strip, (J%(k x )), with N — 20, t\ — 1, 
£2 = 1/3 and M — 0. For clarity, the size of the dots is 
proportional to the fourth power of (Jn(k x )} and the blue 
(red) dots indicate negative (positive) values, (a) p = 7r/3 
showing counter-propagating contributions to the currents. 
The net currents vanish in the bulk but are non-zero close 
to the edges; see Fig. 4(a). (b) p — it /2 showing balanced 
contributions throughout the strip, leading to (J*) = 0. 




Figure 11. Orbital magnetization Ad with N = 20, £i = 1, 
£2 = 1/3 and p = 7 t/3 . Numerically we observe that Ad 
vanishes within the topological phases when M = zbsinc^. 
We also observe that Ad has extrema for M — 0 and 
M — d=\/3sin(^; the latter correspond to the boundaries of 
the topological phases as indicated by the dotted lines. 


Jn{1 ) 



Figure 12. Time dependence of the edge current J^it) after 
a quantum quench from the topological phase with M — 1.4 
to the non-topological phase with M = 1.6, keeping t\ — 1, 
£2 = 1/3 and p = 7r/3 fixed. This corresponds to the main 
panel of Fig. 5 over a longer time duration. The three panels 
from the top to the bottom correspond to strips of width 
N — 20, 30 and 40 respectively. The appearance of resurgent 
oscillations is evident in all three panels due to the finite width 
of our system. The dashed lines indicate the time-scale £ = 
(iV/2)\/3/2c at which signals propagating from the two edges 
meet, corresponding to the onset of finite-size effects. The 
horizontal line corresponds to the ground-state expectation 
value of the edge current for the post quench Hamiltonian. 


Figure 10. Equilibrium edge current J^ for N = 20, £1 = 1, 
£2 = 1/3 and M — 0 with a 5% particle (solid) or hole 
(dashed) doping. In contrast to the half-filled case shown in 
Fig. 4(b), the currents no longer vanish at p = 7t/2. Instead, 
the currents have the same periodicity as the Hamiltonian. 
The increase or decrease of the edge current at p = 7t/2 re¬ 
flects both the sign of the doping and the Chern index. 

new values that are comparable to the ground expec¬ 
tation values evaluated for the final Hamiltonian. In 
these examples we necessarily see finite-size effects due 
to the finite width of the strip. At late times, we see 
resurgent oscillations due to the light-cone propagation 
of currents into the interior of the sample; see Fig. 12. 
The onset timescale for these resurgent oscillations in¬ 
creases with the width of the strip. This timescale is 


given by £ = d/2c where d = N\/ 3/2 is the width of 
the sample and c = 3t\/2h is the effective speed of light. 
This timescale is indicated by the dashed lines in Fig. 12. 
In order to avoid finite-size effects in our predictions we 
therefore restrict the domain of our simulations to be 
within this time interval. The agreement between the re¬ 
sults for N = 30 and N = 40 in Fig. 5 of the main text 
highlights that we are probing the intrinsic dynamics of 
the edge currents, before finite-size effects play a role. 

For completeness, in Fig. 13 we show quenches to the 
topological phase. Numerically we observe that the edge 
currents approach new values that are very close to those 
evaluated in the ground state of the final Hamiltonian. 
Likewise, the oscillation frequencies coincide for quenches 
to the same final Hamiltonian. 
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Figure 13. Edge current J%(t) following a quantum quench 
within the topological phase (filled circles) and from the non- 
topological phase to the topological phase (empty circles). 
We set N — 30, t\ — 1, £2 = 1/3 and Lp — 7r/3 and consider 
quenches of the mass parameter M = 0.5 —>>1.4 (full circles) 
and M = 2.2 —> 1.4 (empty circles). The horizontal line 
corresponds to the ground-state expectation value of the edge 
current for the post quench Hamiltonian. The coincidence 
between the oscillation frequencies is consistent with the fact 
that we quench to the same final Hamiltonian. 







